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In the framework of a single scalar field quintom model with higher derivative, we construct in this 
paper a dark energy model of which the equation of state (EOS) w crosses over the cosmological 
constant boundary. Interestingly during the evolution of the universe w < — 1 happens just for 
a period of time with a distinguished feature that w starts with a value above ~1, transits into 
w < —1, then comes back to m > — 1. This avoids the Big Rip jeopardy induced by m < — 1. 



PACS number(s): 98.80.-k, 95.36.+X 

I. INTRODUCTION 

In 1998 the analysis of the redshift-distance relation of Type la supernova (SNIa) established that our universe is 
currently accelerating [l], Q , which has been further confirmed by recent observations of SNIa at high confidence level 
0, 0, 0] ■ Combined with other observations and experiments this strongly indicates that our universe is dominated 
by a component with a negative pressure, dubbed dark energy, which has been studied widely in the recent years. 
Various candidates for dark energy, such as cosmological constant, quintessence @, phantom [8] and the model of 
k-essence with non-canonical kinetic term [9I [lOj. have been proposed. 

Though the recent analysis on the data from the Supernova, cosmic microwave background (CMB) and large scale 
structure (LSS) show that the cosmological constant fits well to the data, the dynamical models of dark energy are 
generally not excluded, see Ref. for a recent review on dynamical dark energy and actually a class of models, the 
quintom, with an EOS w larg er than —1 in the past, less than —1 today, and evolving across —1 in the intermediate 
redshift is mildly favored [12, • 

There have been many increas ing activities in the studies on the quintom-like models of dark energy 
E, El H iH, HI m, Hi, Hi, Hi, Ha, Hi, Hi, H iH, [slJMJslJjl in the recent years. For the conventional single 



scalar field, the EOS w cannot cross over w = —I 33, [S^, |3^ due to the instabilities of dark energy perturbation 
during evolution. In this paper we focus on a particular example of the quintom dark energy with a single scalar field 
and study its cosmological evolution. 

The model under investigation is a generalization of the model in Ref. [40] with the Lagrangian given by 

C = £(0, X, □0n0, V^Vi.(?!)V^ Vf/)) . (1) 

In the lagrangian ([T]), X = iV^(/)V^0 and □ = V^V^. From ([T]), it is straightforward to get the equation of motion 

1^ - V,(§ V^^) + □(^) + V.V,(f V^V^) = , (2) 
and the energy-momentum tensor 



,dC „ ..„., , _ ,dC_. , „ ,dC. 



Vp( — V^V''(/.)V> - Vp( — VV^) W + Vp{—V^V'cj))VP(b , (3) 



where s is defined as 



1 



(4) 
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This model for the moment is just an effective theory and we assume that the operators associated with the higher 
derivatives can be derived from some fundamental theories, for instance due to the quantum corrections or the non- 
local physics in the string theory 4l[ 4^ [43l|. In addition, with the higher derivative terms to the Einstein gravity, the 



theory is shown to be renormalizable [4^ which has attracted many attentions. Recently, higher derivative operators 
have been considered to stabilize the linear fluctuations in the scenario of "ghost condensation" [4^ . 

In Ref. [13], the authors studied in detail a specifical model which is equivalent to the uncoupled two-field model 
(one is the quintessence-like and the other is the phantom-like) and gives rise to the fate of the universe dominated by 
the phantom-like field, as was implied in two-field quintom model J^2, 14, 15, 16. ITlllsj. Thus the universe might face 
someproblems such as ending in "big rip" (There are some other possibilities of the fate of phantom universe, 
see [13, \^ [131 )• In this paper, however, we generalize this kind of model and show some new possibilities of the 
cosmological evolution. Specifically we will show that in our case there is a possibility that w; < — 1 happens just for 
a period of time. Hence in this scenario it is free of the big rip jeopardy. 

II. MODELS WITH HIGHER DERIVATIVE 

Consider the Lagrangian with the following form 

C = iA(^)V,,^W + §^(°<^)' - ^(<^) ■ (5) 

As proved explicitly in 40], such a model is classically equivalent to a two-field model. This is consistent with the 
general result stated by Ostrogradsky theorem (a nice review on this theorem can be found in 51]), higher derivative 
theory brings more degrees of freedom. For the purpose of seeing how the model in ([5]) is equivalent to two-field 
model, we introduce an auxiliary field x and a new Lagrangian, 



1 M^Y^ 
C = -A(0)V^0V^0 - V,,0V^x - - V{^) . (6) 



Through the variation on x, one can see that 



x = Sfa*. (7) 



By substituting it into ([5]) and dropping out a total derivative, the new Lagrangian ^ reduces exactly to the old 
one ([S]). But the former has the form of two-field without higher derivatives. This is more clearly by further field 
transformations. Consider the field (/> as a function of two independent scalar fields x and ip. Using the transformation 

= (8) 

Eq. ^ can be written as 

C - lAm^ry.^pw^p - iA(0)(|^)2v^xV^x - - vm^p. x)] , (9) 



as long as the transformation ([5]) satisfies 



So the model in ^ is equivalent to the two-field model ([9]). The kinetic terms of ip and x have to take opposite signs 
and one of them must be ghost, which depends on the sign of A{(j)). Furthermore, we can choose </> as a function of 
ip + X and set 

(11, 

So, we get a more elegant form for the Lagrangian 

1 1 M^iX^ 
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The form of (/> as a function of x tp is determined by the conditions (fTO|) and ((TT|| . except constants. However, 
there are some constraints on and C(0). As indicated in (UHl) and pT|) . A{(t)) and C((/)) should not be zero 
in the configuration space. Once the functions A{(f)), C{(j)) and V{(t)) and the initial conditions are determined, the 
evolution of the system in ([5]) is fixed. Based on the discussions above, we can simplify the analysis on the single field 
model (O by analyzing the two-field model (fT2)) . 



III. A SPECIFIC CASE WITH TEMPORARY PHANTOM PHASE 



We will study the evolving behaviors of such type of models for a simple case with A{(f>) = —1, V{cj)) — 0. Thus 
the Lagrangian is simplified as 



= -;;V,*V»,i + -VjxV'x - „^,J' . „ . (13) 



From Eqs. (fTO|) and (fTT|) . one has the relation (p = —{ip + x)- The equations of motion for the two fields are: 



,2 



2C2(V^ + X) 

where the prime is the derivative with respect to ip + x- In the Friedmann- Robertson- Walker universe, □ — /dt^ + 
md/dt with H being the Hubble expansion rate. Furthermore, one can easily get the density and the pressure of 
dark energy, which are: 

„;.2 V.2 A/r2^2 

(15) 

(16) 



and the equation of state is: 



p = 




4 


, M^^,x' 


2" 




2C(V + X) ' 


p = 


2" 


4 


2C(V' + X) ' 


p _ 


1 






p 
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^ 2 2C(V+X) 



(17) 



where the dot represents the derivative with respect to time, x is quintessence-like and ip is phantom-like. They 
couple to each other through the effective potential, i.e., the last term in the right hand of equation P^ : 

M^iX^ 

The coupling function considered in this paper is 

Cicj^iiP + x)] -Co[^ + arctan(^)] = Co[^ + arctan(^^^±^)] , (19) 

thus 



Co A 

Mpill + ^i^ + xf] 



C'{^ + X)= . ,° , ■ (20) 



Because C{ip! + x) is almost constant when |-0 -|- xj 0, V' and x are nearly decoupled at these regimes. We call 
them as "weak coupling" regimes. By contrast, the two fields couple tightly in the "strong coupling" regime where 
\ip + x\ ^ 0- In the weak coupling regime, as shown in Eq. (fT4| . the phantom-like field ip behaves as a massless scalar 
field and its energy density — (l/2)-0^ dilutes as a^^, where a is the scale factor of the universe. The quintessence-like 
field X has a mass term with — Mpi/ y^C{ip + x)- Its behavior is determined by the ratio of m^/H. If <C H, 
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X is slow-rolling and it behaves like a cosmological constant. On the other hand in cases ^ H, the kinetic term 
and potential oscillate coherently and assembly evolve as a^^ , just like that of coUisionless dust. 

It is however more complicated in the strong coupling regime. If the effective potential Veff is negligible in 
comparison with the energy density of the dominant component in the universe, as in the period of radiation and 
matter domination, the Hubble damping terms ZHij: and 3i?x dominates in the equations of motions, see Eq. (jl4p . 
both if] and x ^-re "frozen" and the equation of state of dark energy w ~ — 1. If Veff contributes significantly to 
the energy density of the universe when dark energy is in the strong coupling regime, both of the fields will relax to 
their equilibrium points with large velocities. Because C{ip + x) ^-^d C"(V' + x) ^-^e positive, the acceleration of '0, 

^ca^(^+tf ^^""Ser than that of x- ^2C^HilZd ~ c{l+x) ^' Provided x is positive. Hence, the ratio of ^/-Vx^ 
will increase and from Eq. (fT7|) the equation of state w will become less than —1 when ip'^ > x^. Otherwise, if x is 
negative, w lid "^iH decrease with the expansion of the universe. 

It is necessary to address more about the equilibrium points of V' and x- Through simple algebraic calculations of 
dVeff/dip ~ and dVeff/dx — 0, we find the equilibrium point of ij: is in infinity except for the line of x = 0. In 
the X direction, the zero point x = is the unique minimum. So, if the initial condition is chosen as x is not close to 
zero, the phantom-like field ij) will roll up the potential monotonically, this is along the direction of ^ —oo and to 
reduce C(V' + (/>). The quintessence-like field x will relax down to the zero point. In a word, the system will evolve 
to the weak couphng regime and the field x will get a large mass = MpJ ^JCi^j} -\- x) because C('(/' + x) ^ as 
■0 -I- X ^ — oo. In Fig.l, the effective potential (the green line), x (the black line), and ip (the red line) as functions 
of Ina are plotted. 

We plot the EOS of such a dark energy model in Fig. 2 and Fig. 3 with different initial values. We choose such 
initial conditions that the dark energy starts in the strong coupling regime and x > at high red-shift. We can see 
that the dark energy is frozen quickly at the initial time. Then it becomes significant around the redshift of z ^ 1, 
and the phantom kinetic term ■0^ becomes larger than the quintessence kinetic term x^ . The dark energy crosses the 
boundary of cosmological constant and its equation of state becomes less than —1. In the future, the system will 
evolve to the weak coupling regime, tp^ will dilute quickly and x will get a large mass. The whole system will behave 
as cold matter and its equation of state will oscillate around the point of w = 0. So, the universe will exit from the 
phantom phase and end in the state of matter-domination. This is consistent with the analysis above. Our result 
resembles the late time behavior of "B-infiation" model [s^ . 
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FIG. 1: the effective potential Ve// x Co (the green line), x (the black line), and ip (the red line) as functions of Ina for 
Co/Mli = 6.0 X lO^^^M-^ A = 25, and the initial values are ipi = -0.26Mpi, 7/>i = 3.52 x 10~^^M|;, Xi = 0.25Mp;, Xi = 
-3.62 X IQ-^'^Mli, with noEo- 0.73. 

The cosmological meanings of such a model are as follows: 

• In this model, we give a kind of cosmological scenario which differs from the former quintom dark energy models. 
This model can give twice crossing of and can give a nice exit of phantom which stays just for a period of time, 
avoiding "big rip" to happen. 
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FIG. 2: The equation of state w as a function of Ina for C/M^i = 6.0 x lO^^'^M^^, A — 25, and the initial values are 
1/-, = -0.26Afp,, = 3.52 X IQ-'^'A/^;, x^ = 0.25Mp,, x» = -3.62 x lO'^'^Mli, with Qdeo ~ 0.73. 
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FIG. 3: The equation of state w as a function of Ina for C/Mpi = 6.0 x lO^^^M^J]^, A — 25, and the initial values are 
i>, = -0.26Mp,, = -2.84 X lO-'^^M^,, x» = 0.25Mp,, x^ = 2.74 x lO-'^'M^,, with Qdeo ~ 0.73. 

• This kind of w gives an example of feature EOS in model building in field theory, which is being paid more and 
more attention on by current data- fitting (ssj . 

• At late time, w oscillates with high frequency around zero, so the universe will end accelerated expanding at some 
time. 



IV. SUMMARY 

In summary, we pointed out that for quintom model of single scalar field with higher derivative, the EOS w can 
have some interesting behaviors. We provide a scenario where the EOS w starts from above —1, transits to below 
— 1, then comes back to above —1. This in some sense implies that w < —1 may happen just for a period of time. 
With this behavior of w, this scenario might be able to get rid of the problem caused by phantom. For phantom 
theory, the EOS w is always less than —1, and there might be many problems such as ending in "big rip" as well as 
other singularities, and quantum instability. However in some class scenario of quintom with higher derivative, the 
evolution of iz; < — 1 is very much like the "tachyon" existing only during the phase transition [54| . Thus the "big 
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rip" problem [4^ is solved in such a context. 

The quantum instability of phantom universe is discussed in detail in Ref. [ssl [soj . Our quintom model based on 
higher derivative theory also has ghost degree and the problem of quantum instability. However, as pointed out in 
Ref. [57*1, the problem of quantum instability arises because 4> s-nd '^4' S'l'e quantized in canonical way independently. 
In fact, both of them are determined by cj), and a more appropriate quantization method seems to be possible to avoid 
the instability. 

This work not only provides some new examples for the list of field theory model of dark energy, which deserves 
being included in a more thorough classification, but useful for stimulating the studies on the field theory with higher 
derivatives. 
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